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ABSTRACT 

Self-consistent factorization of two-body residual interaction is proposed 
for arbitrary density- and current- dependent energy functionals. Following 
this procedure, a separable RPA (SRPA) method is constructed. SRPA 
dramatically simplifies the calculations and demonstrates quick convergence 
to exact results. The method is tested for SkM* forces. 
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1 Introduction 



Effective nucleon- nucleoli interactions (Skyrme, Gogny, ...) are widely used 
for description of nuclear properties (see, for example, Ref. fT]). However, 
their application to nuclear dynamics is rather limited even in the linear 
regime. The latter is usually treated within random-phase-approximation 
(RPA) and assumes diagonalization of matrices with the rank determined by 
the size of the particle-hole (ph) configuration space. In deformed and heavy 
spherical nuclei this space is impressive and condemns to a huge computa- 
tional effort. As a result, the effective forces are mainly applied to study 
ground state properties while RPA calculations are very limited. 

RPA problem becomes much simpler if the residual two-body interaction 
is factorized (reduced to a separable form): 



This allows to avoid dealing with high-rank matrices. The main trouble 
is to accomplish the factorization self-consistently, with minimal number of 
separable terms and high accuracy. 

Factorization of the residual interaction is widely used in nuclear the- 
ory but mainly within trivial schemes exploiting one separable term with an 
intuitive guess for the separable one-body operator X. The strength con- 
stant k is usually fitted so as to reproduce available experimental data (see 
e.g. 0). The separable interaction thus obtained does not depend on nu- 
clear density and is not consistent with the mean field. Obviously, accuracy 
and predictability of such schemes are rather low. Several self-consistent 
schemes 3-8 proposed during last decades signified a certain progress in this 
problem. However, these schemes were not sufficiently general. Some of them 
were limited to analytic or simple numerical estimates 3-6 the others were not 
fully self-consistent |7j or covered only particular effective forces [§]. 

In the present paper we propose a general self-consistent separable RPA 
(SRPA) approach relevant to arbitrary density- and current-dependent func- 
tionals. The self-consistent scheme of Ref. |)] is generalized to the case of 
several separable operators. The operators have maxima at different slices 



mnij 



^2 < mn\V res \ij > a^a+ajdi 




X k = ^2< p\X k \h > a+a h . 

ph 
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of the nucleus, which is crucial for accurate reproduction of V res and high 
numerical accuracy. SRPA is tested for the case of SkM* functional |J . 

2 General formulation of SRPA 

The nucleus is assumed to undergo small-amplitude harmonic vibrations 
around HF or HFB ground state. The starting point is a general time- 
dependent energy functional E(Jf(f, t)) depending on a set of arbitrary neu- 
tron and proton densities and currents J"(r, t) (s = n,p; a labels densities 
and currents) 

E(J«(r,t)) =< V(t)\H\V(t) >= J H(r)dr (2) 

where \^(t) > is the wave function of the vibrating system described as the 
time-dependent Slater determinant. Time- dependent densities and currents 
are determined through the corresponding operators as 

J?(r,t)=<*(t)\J:(r)\*(t)>. (3) 

The wave function \^(t) > is obtained from the static HF Slater determinant 
\^o> by the scaling transformation || 

K 

!*(*)>= II exp[-i(q sk (t)-<q sk >)P sk )exp[-ip sk (t)Q sk )\ty > (4) 

k=l 

where Q sk and P sk are, respectively, T-even and T-odd generators fulfilling 
the equations 

Q~sk = Qsk, TQ sk T = Qsk, [H , Qsk] = iPsk, 

P s \ — Psk, TP sk T 1 = —P sk , [H,P s k] = —iQsk (5) 

and q sk {t) and p sk {t) are the corresponding T-even and T-odd harmonic de- 
formations given by 

q sk (t) =<tt(t)|Q afc |*(t)>, P sk(t) =<tt(t)|Afc|tt(t)>, 

<q sk >=<^o\Qsk\^o> ■ 

(6) 
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The HF single-particle Hamiltonian is 

a,s s 

Using the scaling (f|), it is straightforward to find time-dependent density 
and current variations 

6j?(?,t) =< y(t)\j:\^(t) >-< ^ \j:\^ >= 

= £ {(qsk(t)-<q sk >) < \[P S kJ:(m > 

k s=n,p 

+ Psk (t) < \ [Q sk J«(r)}\ >} (8) 
and the response Hamiltonian (vibrating single-particle potential) 

h(f,t) ~h (r)+h res (r,t), 

h res (r,t) = B^#]*# (F,f) =EE[^1^(^)*W 

a's> OJ s' as a ' s ' 0J s 0J s' 

= Y,{(Vsk(t)-<q s t>)X sk (r)+p sk (t)Y sk (f)}. (9) 

sk 

In the above equations we introduced one-body operators 

s'a'a s' s 

where X = X + is T-even, Y" = Y + is T-odd, and the notation 

<[A,B]>=<* \[A,B]\* > (11) 

is used. The moments of the operators are 

< 5X sk (t) > =< y(t)\X sk \V(t) >-< ^ \X sk \^ >= 

= iJ2(^'k'( t )- < Qs' k '>) < \[P s ' k ',X sk ]\ >= 

s'k> 

= £(9«w(*)-<^k'>)«7fcW (12) 

s'fc' 
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< 8Y sk (t) > =< *(t)\Y ah \V(t) >-< V \Y sk \V >= 
= iJ2p s 'k'(t) < \[P s ' k ',Y sk \\ >= 

s'k' 

= ^2Ps'k'{t)v7k)s'k' (13) 
s'k' 

where we introduced the inverse strength matrices 

K sk,s'k' = K s'l',sk = ~i< \[X,k, Ps'k'W >= (14) 

= / < \[JT,p*\\ > < \[%,P**]\ >, 

OLd 1 s s 

V7k,s'k' = V7'l',sk = -i< \[Ysk, Qs'k'W > (15) 

= / ^£[^^] < \[J?,Q.k]\s > < P^Qs'k'}\ > ■ 

act' s s 

Harmonic collective shifts q s k(t)—<q s k> and velocities p sk (t) read as 
q s k(t)-<q sk >= q sk cos(ujt) = ^q sk (e iujt + e~ iujt ), 

p sk (t) = p sk sin{ut) = yPakie™ - e^*). (16) 

Following the Thouless theorem, the perturbed wave function (0) can be 
also written as 

\^(t) >~ (1 + J2eph(t)ala h )\y > (17) 

ph 

where c p hit) are the time-dependent particle- hole contributions to the given 
excited state. Substituting (|P7|) into the time-dependent HF equation 

ij t \V(t) >= (h + h res (t))\V(t) > (18) 

and using 

Cphit) = c p V^ + c- h e~^\ (19) 
one gets the relation between ci and collective deformations q sk and p sk 

± 1 J2 k [q s k < p\X sk \h > Tipsk < p\Y sk \h >] 

Cphts ~ 2 7^±7j ' (20) 



where e p h = e p — Eh } £ P and Eh are the particle and hole energies, respectively. 

In addition to Eqs. (0)"dH)' the moments < 5X sk {t) > and < 5Y sk (t) > 
can be also expressed in terms of the Thouless wave function ([16|): 



<5X sk (t)> =J2(c ph (t)* <p\X sk \h> +c ph (t) <h\X sk \p>), 

ph 

<5Y sk {t)> =J2(c P h(t)* <p\Y sk \h> +c ph (t) <h\Y sk \p>). (21) 

ph 

Then, using Eq. (^) an equating the moments (|T2"D-(|I3"D and fl2"T|), we obtain 
the system of equations for unknowns q sk and p sk : 

^2Qs'k'[S s 'k',sk(XX) — Kg'lt^] +'^2p s 'k'S s 'k',sk(XY) = 0, 
s'k' s'k' 

J2^ , k' S s'k',sk(XY)+J2Ps'k'[Ss'k',sk{ YY )~V7'k',sk} = ( 22 ) 



with 



S s r k > >sk (XX) = -a— — j{< p\X s > k >\h >*< p\X sk \h > (s ph - u) 
ph £ P h ~ ^ 

+ < p\X s i k >\h >< p\X sk \h >* (e ph + a;)}, 

S s ' k ' !Sk (YY) = Y — P\ Y s'k'\ >*< P\ Y sk\h > (e ph - u) 
P h £ P h ~ ^ 

+ < p\Y s ' k '\h >< p\Y sk \h >* (e ph + uj)}, 

S s , k , tSk (XY) =~i^2— o{<p\X s ' k '\ >*<p\Y sk \h> (e ph -u) 

ph £ P h ~ u 

+ <p\X s , k ,\h ><p\Y sk \h >* (e ph + u)}. (23) 
The matrix of the system (p2|) is symmetric and real. By equating the deter- 



minant of the system to zero, detD = 0, we get the dispersion equation for 
eigenvalues uj v . 

It can be shown that the same equations can be obtained through the 
standard RPA equations with the Hamiltonian 

Hrpa = h + V r es, (24) 

where 

Vres = — - X! Y,[^sk,s'k'X sk X s / k i + fj sk}S i k iY sk Y s i k i}. (25) 
sk s'k' 
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The strength constants k and 77 are obtained by inversion of the matrices 
© and (H). 

Eqs. (0), (0), (0), 
tions. 

The rank of SRPA matrix ( p2[) is 4.K" with K = 3 — 6 (see next section). 
So, high-rank matrix of exact RPA shrinks to low-rank SRPA matrix. This 
minimizes the computational effort. At the same time, as is demonstrated 
below, SRPA provides high accuracy of the calculations. 

The approach is fully self-consistent and does not need any adjusting 
parameters in addition to the starting functional. Unlike the trivial separable 
schemes, SRPA gives analytical expressions for both separable operators and 
strength constants. Relative contributions of different separable operators 
are self-consistently regulated for every RPA state. 



20), (22), (23) constitute the set of SRPA equa- 



3 SRPA test with Skyrme functional 

SRPA was tested for the particular case of SkM* functional |§. Isoscalar 
and isovector densities (nucleon, kinetic, spin and spin-orbital) and currents 
have been involved. The strength function 

b L (XX,gr^u) = < u\M(XX)\0 > 2 p(u - u v ) (26) 

with Lorentz weight p(x — y) = ^ ( x ~- y ) 2 -(A/2) 2 (A = 1 MeV is the averaging 
parameter) has been calculated for isoscalar E2 and isovector El transitions 
in 40 C a and 208 Pb and compared with the results of exact RPA. In (|26|). 
M.(X\) is operator of XA-transition. 
The set of K = 2 generators 

Q k (r) = R k (r)(Y^ + h.c), P k = i[H, Q k ] (27) 

has been used with Ri(r) = r A , i?2( r ) — r A+2 . The first (^-generator 
coincides with the corresponding i?A-transition operator in the long-wave 
approximation, and the other creates operators X sk and Y sk with maxima in 
the interior. 

The comparison between SRPA and exact RPA results is exhibited in 
Fig. 1. Already first operator gives excellent description of E2 giant reso- 
nance. However, it is not enough to reproduce El data. The second operator 
considerably improves the agreement and provides a satisfactory description. 
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4 Conclusions 



A general procedure for self-consistent factorization of the residual interac- 
tion is proposed for arbitrary density- and current- dependent functionals. 
The separable RPA (SRPA) constructed in the framework of this approach 
dramatically simplifies the calculations while keeping high accuracy of nu- 
merical results. The latter was demonstrated for the case of SkM* functional. 

SRPA can be used for description of EX and MX response in both spher- 
ical and deformed nuclei. In the case of deformed nuclei it allows to take into 
account the coupling between electric and magnetic modes. The approach 
can serve as a good basis for anharmonic corrections, description of vibra- 
tional states in odd and odd-odd nuclei. SRPA can be especially useful for 
investigation of dynamics of deformed nuclei where applications of effective 
forces (Skyrme, Gogny, ...) are very scarce. One of the most promising lines 
of future studies is dynamics of exotic nuclei obtained in radioactive beams. 
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FIGURE CAPTIONS: 



Figure 1. Isoscalar E2 and isovector El strength functions in Ca and 
208 Pb. The lines represent exact RPA results (full) and SRPA ones with one 
(dotted) and two (dashed) operators. 
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